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ABSTRACT 



Orifice and louver sound power spectra are investigated, using an 
intensity probe, at various wind speeds in a low noise, semi-anechoic , 
subsonic wind tunnel for free stream velocities below 50 meters per 
second. The radiated noise is created by turbulent flow over various 
orifice and louver geometries which are flushed mounted into the wall 
of a long duct. Five orifice samples of rectangular shape and various 
transverse dimensions as well as four louver samples with multiple 
rectangular and circular orifices are tested. Also investigated is 
the effect of the leading and trailing edge angle on the radiated 
sound power. The scaling laws of the excitation frequencies and the 
speed/power laws are presented for ratios of the boundary layer thick- 
ness to the transverse orifice dimension from 1.01 to 4.29. A 
detailed theoretical model is developed for rectangular shaped aper- 
ture orifices and louvers based on the work by Ffowcs Williams [9], 
Nelson [11], and Corcos [12]. The model describes the predominant 
variables which effect the radiated flow noise. To ascertain the val- 
idity of the model, it is used to collapse the experimental power 
spectra for two orifice and one louver geometry. The model showed 
excellent agreement with the actual sound power measurements. Direc- 
tivity plots are also presented to further describe the orifice/louver 
source . 
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1. INTRODUCTION 



1.1 Background 

Turbulent flow over an orifice or louvered openings can generate 
high levels of broadband noise. This phenomenon has been known for 
some time [1] and early investigators like Dunham [2], Rossiter [3] 
and East [4] attempted to identify and control the high vibrational 
levels and possibly the radiated sound levels caused by the interac- 
tion of the flow over these openings. Recent work by Elder, Farabee 
and DeMetz [5] concentrated on the area of cavity resonance where the 
opening provides the feedback path between the closed acoustic cavity 
and the source, the turbulent flow. When the typical dimension of the 
cavity is large compared to the acoustic wavelength, the source can 
excite various room modes of the cavity. If the turbulent vortex 
shedding frequency coincides with one of these room modes, then a 
strong feedback will occur, resulting in an acoustic resonance of that 
mode . 

If the cavity's volume is further increased, multiple room modes 
will be excited and the strong feedback mechanism will be lost. This 
results in a much broader acoustic response and at lower amplitude 
levels than the resonant case which is typically characterized by 
large amplitudes at very discrete tonal frequencies. In practice, the 
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detection of resonating cavities is often easy due to their strong 
discrete tones. However, once these resonances have been corrected, 
or presumably prevented in the design, identifying the cause of the 
remaining flow noise becomes much more difficult. Possible contribu- 
tors to the broadband flow noise are: (1) the turbulent boundary layer 
direct transfer of energy to the structure of the enclosed acoustic 
space, (2) the power flow through the opening (open window effect) 
caused by the acoustic field created by the upstream and downstream 
turbulent sources in the boundary layer, (3) the power flow into the 
acoustic space caused by the source created in the opening due to the 
driving aerodynamic field and the shear- layer edge interaction with 
the downstream edge of the orifice and (4) the vibrational excitation 
of the internal structure caused by the orifice source. The amount of 
noise which the latter contributes depends on the details of the 
internal structure and the orifice/louver attachment design. There- 
fore to predict the contribution from internal structure vibration, 
one must first determine the orifice source. 

1.2 Content 

The contribution of the first three possible contributors to the 
total power radiated into a closed acoustic space is analyzed based on 
the results from several experiments. It is shown that indeed the 
orifice source provides both the spectral character and the predom- 
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inant amount of the total sound power. By analyzing and comparing the» 
results with a derived equivalent source model, the nature of the 
source and the parameters which affect its strength and spectral char- 
acter are identified. The analysis comes directly from experimental 
measurements taken in a subsonic wind tunnel at the M.I.T. Acoustic 



and Vibration Laboratory. 
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2. SOUND INTENSITY MEASUREMENTS 



2.1 Why Measure Sound Intensity? 

To describe an acoustic source, one must measure a quantity which 
is not dependent on the environment or the surroundings. The quantity 
most often identified to describe a source is its power. Knowing the 
source's power, its interaction in any known environment or any known 
surrounding can be determined. 

Until recently, most all experimental measurements to describe an 
acoustic source were accomplished by measuring sound pressure. This 
quantity, however, was the result of a source's interaction with a 
specific environment and utilizing the experimental results in other 
applications was often difficult. One technique that was developed to 
overcome the difficulty in relating one environment to other was to 
isolate the source of interest and to simulate a free field environ- 
ment in an anechoic chamber. Here the radiated sound pressure in air 
was measured, then related directly to the sound intensity. This 
technique has as its biggest limitation, the availability of an effec- 
tive anechoic chamber, particularly for use in low frequency analysis 
at the low frequencies. 



Another technique is to measure the intensity directly utilizing a 
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sound intensity probe. The theory of intensity measurements comes 
directly from two fundamental properties of acoustics. Those are the 
conseirvation of momentum and the definition of intensity. Conserve* 
tion on momentum states that: 

pdv^/dt - - dp/dr (2.1) 

which can be solved for the radial component of the particle velo- 
city, Vj- 

Vj- - . (1/p) J dp/dt dt (2.1a) 

A finite difference approximation is used to calculate the pres- 
sure gradient using two identical microphones by 

dp/dr - lim (Ap/Ar) a (p^-pg)/Ar (2.2) 

Ar-^0 

where A and B correspond to channels A and B for the respective 
microphone inputs and Ar is the finite microphone spacing. The 
instantaneous pressure is taken as the mean value of the two pressure 
signals , 

P “ (Pa+PbV 2. (2.2a) 



The finite difference approximation is valid provided the separa- 
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tion distance, Ar, is small compared to the acoustic wavelength. 

Thus, a finite difference approximation error becomes significant at 
high frequencies. For a plane sinusoidal wave, propagating along the 
axis joining the two microphones, the finite difference approximation 
assumes the free field phase between the two microphone positions to 
be equal to kAr, whereas the actual measured free field phase is 
Sin(kAr). Therefore, the measured intensity, I', is related to the 
actual intensity, I, by 

r/1 - Sin(kAr)AAr (2.3) 

and the resulting measured Intensity underestimates the actual 
intensity by at high frequency by 

Lg - 10 log(Sin(kAr)/kAr) . (2.3a) 

All the intensity measurements in this thesis were obtained with a 
spacing of Ar - 50mm. Thus the finite difference error for this 
spacing, at the upper analysis bandwidth frequency of f - 1600 Hz, was 
Lg - - 1.7 dB, 

From the definition of intensity and by substituting equations 2.2 
and 2.2a, the measured sound intensity becomes 



I 
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Ir " PVr = - ( (PA+PB)/2Mr)/(pys^-pB)dt. (2.3b) 

By measuring the intensity over a control surface which surrounds 
the source, the power can be calculated from: 



P - / Ir dA (2.4) 

A 

The advantages of measuring intensity directly are: (1) the 

effects from diffuse noise contaminating the experimental results are 

eliminated since the integral of this noise intensity around any 

closed control surface is zero, (2) the sound power calculated 

describes the acoustic source in a form which is readily usable to 

design engineers, (3) measurements can be made in the source's near 

field and (4) there are no restrictions on the shape and size of the 

control surface. A simple yet effective analogy to describe the 

rationale in measuring sound intensity, thus power directly follows: 

In attempting to design a heating system, an engineer 
solicits specifications from two companies. One com- 
pany describes its heating system by stating the 
electrical power and the thermal efficiency of the 
unit. The other company describes the temperature at 
various locations inside an experimental enclosure at 
steady state conditions. 



It is obvious in this example that describing the temperature 
field within one environment is practically useless to the design 
engineer whereas knowing the power and thermal efficiency, the 
designer can determine the interaction of this source with his partic- 



ular environment. 
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In the experiments which follow, every attempt was made to uti- 
lize the advantages of measuring intensity directly and simulating a 
free field environment by taking the measurements in a semi -anecholc 
chamber, thus minimizing any enclosure interactions which could affect 
the radiated sound power. 

2,2 Ideal Versus Practical Sound Intensity Analyzer 

2.2.1 Active and reactive sound fields 

Sound fields exist in two parts: an active part which transports 
sound energy by virtue of the fact that the sound pressure and the 
particle velocity are in phase and a reactive part which stores sound 
energy and in which the sound pressure and particle velocity are in 
quadrature. A sound field is said to be more reactive the more energy 
is stored relative to that transported. Sound field reactivity is 
defined as the difference between the measured sound pressure level, 
Lp, and the maximum intensity level Lj^. 

2.2.2 Ideal sound intensity analyzer 

An ideal sound intensity analyzer would measure only the active 
part of the total sound field and would Indicate sound intensity val- 
ues approaching 0 watts/m^ in a highly reactive (highly diffuse) 
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sound field. 

2.2.3 Practical sound intensity analyzer 

A practical sound intensity analyzer can never display such a low 
intensity level due to the small but ever present phase-mismatch 
between the two microphone channels. It is this phase mismatch which 
accounts for the lower limit in the dynamic range of a practical sound 
intensity analyzer. 

2.3 Determining the Limits of Error for the Measured Intensity 

Errors are introduced in any practical intensity measurement by: 

(1) the influence of the phase mls-match in the probe, and 
the analyzer, 

(2) the Influence of the phase of the sound field of known reacti- 
vity, and 

(3) using a particular spacer between the two microphones. 

The limits of these errors can be determined using a phase- 
reactivity nomogram [see reference 6). This nomogram graphically 
determines, for a given measured reactivity (Lj^-Lp) and a 
particular spacer, Ar, according to the relationship: 
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Li - (Lp+0.16) - 10Log(1000/Ar x C/f x <^f/360] (2.5) 

The phase of the sound field, together with the phase of the probe 

and the analyzer, ^f ^a ^ ^m* limits of 

error in the measured phase at a particular frequency. From these 
phase limits, the nomogram provides the limits of error in dB for the 
measured intensity. 

2.4 Sound Intensity Probe 

The specific details on the operation of the intensity probe used 
in this experiment are contained in reference (7). Briefly, the 
intensity probe consists to two B&K phase matched, 1/2" condenser 
microphones, type 4177/4165, separated by a spacer 50 mm in length and 
on a common axis. One microphone is a reference while the other 
microphone provides the pressure gradient measurement. 

2.5 Intensity Measurements 
2.5.1 Calculating intensity 

Sound intensity measurements were analyzed using a h&K 2032 Spec- 
trum Analyzer. The sound intensity was calculated according to the 
relationship : 
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I(k) - - (l/«Arp) Im[GAB(k)) 



( 2 . 6 ) 



where Ar is the microphone spacing in meters, A corresponds to 
Channel A which is the reference microphone input, B corresponds to 
Channel B which is the second microphone input, k is the frequency 
index and p is the density of air, calculated according to the rela- 
tionship ; 



where p is the pressure in mbar and T is Che temperature in degrees 
centigrade , 

The cross spectrum, G^(k) , is the complex ensemble average of the 
complex products of the complex conjugated one-sided instantaneous 
spectrum Gy^(k) and the one-sided instantaneous spectrum Gg(k) . 



The instantaneous one-sided spectrum for channel A or B, G^(k) , 
Gg(k) was calculated according to: 



p - 0.348534 p/(273+T) 



(2.7) 



Gab(^) - GA(k).GB(k) 



( 2 . 8 ) 



' $A(k) for k - 0 

G^(k) - - 2 S^(k) for 1 < k < N/2 - 1 



(2.9) 






0 



for N/2 < k < N-1 
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where N is the number of samples per ensemble. For this experiment, 

N was equal to 1024. Sy^(k) is the two-sided instantaneous spectrum 
for channel A which is determined from the recorded time record, a(n) 
according to: 

SgOc) -^(Q(n).a(n)) (2.10) 

where ^ is the forward discrete fourier transform and w(n) is the 
time record weighting function. In this experiment, a Hanning window 
weighting function was utilized. 

The forward discrete fourier transform, ^ , is defined as: 

N-1 

S^(k) -1/NS a(n)*exp(-j 2nkn/N) (2.11) 

n-0 

where n equals the time index times the sampling interval, nAt 
and k equals the frequency index times the frequency resolution, kAf. 
The frequency resolution for this experiment was 2Hz for an analysis 
bandwidth of 1600 Hz. 

2.5.2 Measurement techniques 

Two techniques can be used to measure sound intensity using an 
intensity probe. They are: 



. 
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(1) continuous spatial averaging 

(2) discrete point averaging 

Continuous spatial averaging consists of sweeping the intensity 
probe at a controlled rate over each surface while continuously samp- 
ling and averaging data. This technique produces the most accurate 
spatially averaged results but suffers slightly in repeatability. 
Several repeated runs, often four or more, were required to obtain the 
required repeatability of total sound power (within 0.5 dB on at least 
three runs). Figures 2.1, 2.2 and 2.3 are the results of three such 
repeatability experiments with 0.5 dB or less variation in total sound 
power. Shown is the sound power spectral level from a rectangular 
orifice, 4 1/2" by 1 2/3", and an analysis bandwidth of 2 Hz. 




Figure 2.1: Louver Li sound power spectra repeatability. Shown are four succes- 
sive spectral measurefnents of sound power. 
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Figure 2.2: Orifice 01 sound power spectra repeatability. Shown are three suc- 
cessive spectral measurements of sound power. 




Figure 2.3 Orifice 01 sound power spectra repeatability. 

sive spectral measurements of sound power. 



Shown are five succes- 
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Discrete point averaging consists of dividing your control surface 
into numerous discrete areas and measuring the sound Intensity at the 
center of each discrete area (point) and summing over all areas. This 
technique offers maximum repeatability but at reduced spatially- 
averaged accuracy. The two techniques where both utilized initially 
and compared. The repeatability of the total radiated power using the 
discrete point technique, with 45 points, was 0.3 - 0.5 dB between 
successive runs. These repeatability values were selected as the 
criteria for the continuous spatial averaging technique. The results 
showed that the two techniques gave similar sound power spectrums and 
similar total sound power levels within an accuracy of 0.3 dB. For 
this reason, the continuous spatial averaging technique was selected 
due to its increased accuracy and reduced sampling time. 

The control surface was constructed of 1/4" balsa wood enclosing a 
total surface area of 0.542 square meters. Intensity spectra over the 
frequency range of 52-1600 Hz were taken and recorded on a B&K 2032 
Spectrum Analyzer. Microphone calibration was performed using a pis- 
tonphone and daily temperature and barometer measurements were taken 
to calculate the air density. To ensure that the two microphone were 
indeed phase matched, their phase was measured in a random noise field 
and compared. The microphones were placed at the open end of a long 
cylinder, equal distances from a speaker placed at the opposite end of 
the cylinder. This apparatus was similar in design to a calibration 



pistonphone, type ^200. Figure 2.4 shows the phase difference measured 
over a frequency range of 2560 Hz and clearly demonstrates their com- 
patibility over the frequency band of interest, 52 - 1600 Hz. 




Figure 2.4: Phase difference between two microphones and their pre-anpl if iers 

(channels A and 6) 

In Section 5.2, the errors in dB resulting from these phase mis- 
matches are calculated and shown to be on the same order of magnitude 
as the inherent statistical fluctuations. 
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2.5.3 Power spectrum contamination 

Data below 52 Hz were ignored due to contamination by the 1,0,0 
room mode (38 Hz) of the semi*anechoic chamber. Also at this low 
frequency, the effectiveness of the anechoic treatment was reduced. 
The peak at 462 Hz and its harmonic at 924 Hz which is evident in all 
the spectra was due to the resonance of the cross -duct mode of the 
test section and not at all associated with the orifice source. For- 
tunately, the sound level and narrow bandwidth associated with this 
modal resonance contributes insignificantly to the total sound power 
caused by the orifice source. Frequencies above 1600 Hz contributed 
nothing to the total sound power and the spectral shape followed the 



same trend as from 1200 - 1600 Hz. 
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3 . THEORY 

3.1 Introduction 

The development of a theoretical model to determine the scaling 
laws for orifice/louver flow noise was first accomplished in 1972 by 
Ffowcs Williams [8]. He treated the model problem of flow noise gen- 
erated by a body of turbulence near an infinite, rigid, thin screen 
with circular orifices. His assumptions were: a) the size of the ori- 
fice is much smaller than the distance between them, b) variations in 
pressure within individual apertures are neglected. While these 
assumptions are rather restrictive in their application to many indus- 
trial louvers, the theory does provide a fundamental result: that the 
sound radiated by orifices and louvers can be directly related to the 
mass flow driven through the aperture by the aerodynamic field. 

This fundamental result was used by Nelson [9] to develop a model 
for flow noise generated over practical perforated screens with circu- 
lar orifices. Due to the normal component of the turbulent fluctua- 
tions above the louver, flow is driven through the orifice resulting 
in a fluctuating volume flow through the upper surface. This fluctu- 
ating volume flow constitutes an equivalent monopole source at the 
upper surface. Since the volume flow into the upper surface of the 
orifice creates an equivalent outflow through the lower surface, there 
is an equal and opposite equivalent monopole at the lower surface. 
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The net field which radiates out of the louver and into the closed 
acoustic space is thus the combination of the lower monopole source 
and the contribution from the upper monopole source which is dependent 
on the degree of acoustical transparency of the louver. If the louver 
is acoustically opaque, then only the lower monopole source will radi- 
ate into the acoustic space. If the louver is acoustically transpar- 
ent, then the upper monopole source will radiate into the space, hav- 
ing a cancelling effect on the field of the lower monopole source. 

This would result in a net radiated field which exhibits properties of 
a dipole source. 

Nelson's equivalent source model gave an excellent description of 
the flow noise produced by perforated facing materials used in prac- 
tice even though he assumed a rather simplistic form for the driving 
aerodynamic velocity. His assumptions were: a) the wavelength for the 
sound is much larger than the orifice radius, b) there are many aper- 
tures per acoustic wavelength and no interactions of the flow occurred 
between adjacent apertures which can modify the aperture's conducti- 
vity. In the sections that follow, the theory developed by Nelson and 
Ffowcs Williams will be utilized to model the flow noise created by 
turbulent flow over orifices and louvers. Once the equivalent source 
model has been developed, it will be used to collapse the experimental 
sound power spectra thus providing a verification to the theory's val- 
idity. The coordinated system used for this analysis is shown in 
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Figure 3.1. 




Figure 3.1: Coordinate system for analysis 
3.2 The Equivalent Source Model for Orifices and Louvers 

Ffowcs Williams [8] and Leppington and Levine [10] developed an 
expression for the reflection coefficient, R, for a plane wave of 
single frequency, (j), incident at a polar angle, on an infinite 
louver. For low frequencies and orifice spacings that are large com- 
pared to the orifice size, the reflection coefficient, R, is 

R - l/(l+4aiNAcos^) (3.1) 

where a is the orifice radius, N the number of orifices per 
unit area of louver and k is the acoustic wavenumber, k - o>/Co , 

The field produced by the lower monopole at yj^ can be determined 
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using the principle of reciprocity. The field at yy due to a distant 
source at xj^, is given by p^(l+R) where pj^ is the Incident plane wave 
striking the lower surface of the louver. Thus one can write 

P(yfc) ” (1+R) Q(Xk) exp(iklyj^-xi^| )/4x|yj^-xi^| (3 -2) 

where Q is the strength of the distance monopole at xj^. By 
reciprocity, the field at xj^ due to the monopole at yj^ is given by 
interchanging the source and receiver position vector in equation 3.2. 
Equation 3.2 can be generalized to related the spectral density of the 
acoustic pressure field, Sp(xj^,w) to the spectral density of the 
source Sq(yj^,w) such that 

Sp(xi,,«) - |l+R|2.SQ(Yk,w)/16x2r2 (3.3) 

where r^ - Ixj^-yj^l^, Figure 3.2 is a plot of |r| 2 against the ratio 
(4aN/kcos^) . 




Figure 3.2: The infinite screen reflection coefficient as a function of the ratio 
(4aN/kcos0>. From C9, Nelson, P.A., 1982]. 
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From this plot, It can be seen that, 

|r| 2 - 1 for (4aNAcos$) s 0.3. 

This occurs for open area ratios on the order of 0.01 (see refer- 
ence 11). For all the orifice and louvers tested in this study, the 
open area ratio is less than 0.005 so that 1 r|^ - 1. To verify this 
result for the reflection coefficient, since the assumption of many 
apertures per acoustic wavelength is not valid for the single orifice, 
a separate analysis was performed. Here, in two dimensions, the plane 
wave of single frequency, w, is incident normal on an infinite, rigid 
surface with a single opening, of dimension, 2b. Again the principle 
of reciprocity is used to predict the radiated pressure from the 
single orifice. Using matched asymptotic expansions, for b/A « 1, 
the lowest order solution for the reflection coefficient is again, 

1r| 2 - 1. The details of this analysis are shown in Appendix 1. 

Thus , 

Sp(xi<.,w) - Sq(yi<.,w)A»r2r2 (3,4) 



The instantaneous source strength can be written in terms 

of the weighted integral of the unsteady normal component of velocity 



I 



I 



I 



31 



over the area of the aperture. To take advantage of the homogeneous 
and stationary properties of the velocity in the plane of the wall, we 
can write 

q(yk*t) - - iwpj Jv(yk.t)K(sj^-yij)dA(si^) (3.5) 

00 00 

where yy^ is the vector position of some reference point on the aper- 
ture surface and K(sj^-y^) is the weighting function, i.e., the 
function which converts the integral over the aperture area, sj^ - 
S]^,S 2 » to an integral over all space, defined as 

K(sk*yk) “ 1 ^k'yk inside of the aperture area, A 
- 0 for Sk • yk oustide A 

Following Corcos [12] the spectral density of the source strength, 
SQ(yk»ci>) can be defined by 

JS.^^(zi^-sj^,w)K(zi^)K(si^)dA(zi^)dA(sj^) (3.6) 

00 « 

where Zk is a dummy space vector. Define ^k 
^k " ^k * ^k 

thus 

SQ(yk»^) " S JS^(€k»^)^(Sk)^(sk+^k^^(^k)^(^k) 

00 00 

- p^<j>^ J S.^(«i^,w)5 («j^)dA(fj^) (3.6b) 

00 

where the function, 5(fk)i defined as the orifice response function, 
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is determined wholly by the geometry of the aperture and can be evalu- 
ated as 



“ /K(si^)K(si^+«i^)dA(sjj) —lx area of overlap 

- <x> 

“ (L- I €3 I ) (b- I €i I ) for |€ 3 | < L, |e]^| < b 

- 0 elsewhere 

for the rectangular shaped orifice shown below 




(3.6c) 



Finally, we assume the cross-spectral density, S^(«j^,w), to have 
the product form 

Sw(fk.‘^) - S^(w)A(Kei)B(K€ 3 )ei*^«l (3.7) 

where S^(w) is the spectral density of the normal velocity fluctua- 
tions, K is the convective wavenumber given by K - where is 

the convection speed of the boundary layer eddies and A(Ke^) and 
B(Kc 3 ) are the amplitude functions for the normalized longitudinal and 
lateral cross-spectra, respectively. This form is similar to the 
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expression given by Corcos [12J for the cross spectral density of the 
boundary layer pressure fluctuations. The longitudinal cross - spectra 
of the eddies as they convect in the direction also contains the 

phase term exp(lKc^). A(Ke]^) is plotted in Figure 3.3 as a function 
of its argument for various ratios of separation distance to momentum 
thickness (rj^/5*) for wall pressure. Since the wall pressure is 
related to the volume integral of the normal component of the velo- 
city, wall pressure cross-spectra data can be used to predict the lon- 
gitudinal and the lateral cross-spectra of the normal component of 
particle velocity. Comparing the wall pressure correlation length 



scales from [12] to the velocity correlation length scales from [17], 
we find that both the longitudinal and lateral length scales agree to 
within an order of magnitude, thus supporting the previous statement. 




Figure 3.3: Amplitude function for longitudinal cross -spectra I density of wall 
pressure. From (16, B8N report No. 4110366, 1966] 
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For the various orifice and louver samples tested, at a maximum 
value for equal to 1 , A(Ke]^) varies from a maximum value of 1 to 

a minimum value of 0.6 for orifice 02 at U-20m/s and f-250 Hz (the 
frequency at which the spectral level is lOdB down from the peak 
level). Throughout most of the samples tested, A(Ke]^) - 1 over the 
frequency band of interest. Thus the longitudinal correlation lengths 
were typically several times the length of the orifice aperture. Thus 
the assumption is made that no decay in the eddy structure occurs in 
the longitudinal (y^) direction. 

The normalized lateral cross-spectra in the y 3 direction is given 
by the term B(Kc 3 ) which is assumed unity for no separation and to 
decay very rapidly with increased separation. To justify this assump- 
tion, B(Kc 3 ) is plotted in Figure 3.4 as a function of its argument 
for various ratios of r 3 / 6 *. From this figure, the typical correla- 
tion length, i, is calculated to be on the order of 2 cm. The ratio 
of the lateral correlation length, i, to the aperture length, L, is 
0.18 thus supporting the initial assumption of rapid lateral decay. 

The correlation lengths are determined at the point where the ampli- 
tude function equals 0.4 and at the frequency corresponding to the 
peak power spectral level. 
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Figure 3.4; Amplitude function for lateral cross • spectra I density of wall pressure 
From [16, BBN Report No. 4110366, 1966) 



Substituting equation 3.7 into equation 3.6b and evaluating the 
integral over the rectangular shaped orifice, assuming S^(yj^,o>) to be 
uniform over the orifice area gives: 



S^(yj^,w)*Li»(2/K2)(l-cos(Kb)) (3.8) 

where L is the orifice length in the direction. See Appendix (3) 
for the details of the integration for the rectangular shaped orifices 
and reference [9] for the results for circular shaped orifices. 
Equation 3.8 states that the monopole spectral density is directly 
related to the normal velocity spectral density. 



36 



3.3 Scaling Laws for the Radiated Sound Power 

Having determined the form for the source strength, we can now 
develop the scaling laws for the radiated sound power from various 
rectangular shaped configurations. Assuming that each orifice radi- 
ates sound incoherently, the total source strength for N apertures 
is simply N times the individual source strength for each aperture 
defined by equation 3.8. Substituting equation 3.8 into equation 3.4 
gives : 



The assumption is made that the convection velocity, U^, is equal 
to one half the frees tream velocity, U«. Define the turbulent 
intensity, I, and the Strotihal number, St(i) and St(b) by: 



Finally, integrating equation 3.9 over the bandwidth Ao) with the above 
expressions substituted, yields 




(3.9) 




(3.10) 



St(i) - fi/U«. St(b) - fb/U« 



(3.11) 



5 



(p2(Xk.t))^^ - Np2L.(l-cos(Kb))St(i)(l2)^^U^2fxr2 



(3.12) 
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For spherical waves, the sound power, P, can be determined by; 

5 

(P(t))A« - NpL (l-cos(Kb))St(i)(l2)^ U„/fCo (3.13) 

Equation 3.13 is the final result for the flow noise created by tur- 
bulent flow over orifices and louvers with rectangular shape as 
described in Section (4.2). The controlling parameters for scaling 
sound power are thus: 

(1) Strouhal number based on the lateral correlation length, i 
and aperture transverse dimension, b 

(2) the normal component of the Turbulent Intensity 

(3) the fifth power of the freestream velocity 

(4) the shape factor (l-cos(Kb)) 

In Section (5.10), this model is used to collapse the experimental 
sound power data for orifices 02 and 03 and louver LI. It is clear 
from the collapse that equation 3.13 does indeed identify the predom- 
inant controlling parameters for rectangular shaped orifice and louver 
flow noise. Note, for circular shaped apertures, the above scaling 
factors remain the same except the shape factor is replaced by the 
function, H]^(Kd) where is a Struve function of the first kind and 
d is the aperture diameter. See reference 9 for the details in cal- 
culating equation 3.8 for a circular aperture. 
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4. EXPERIMENTAL APPARATUS 
4.1 Semi-anechoic Wind Tunnel 

The experiments were carried out using the low turbulence, sub- 
sonic wind tunnel of the Acoustics and Vibration Laboratory at MIT. 

The wind tunnel details are shown in Figure (4.1) with the basic con- 
struction characteristics described by Hanson [13] and modifications 
described by Shapiro [14]. The wind tunnel is of open circuit design 
with the test section enclosed in a blockhouse. The tunnel consists 
of a set of flow straighteners , screens, a settling chamber and a 20:1 
contraction leading to a square test section 38cm by 38cm. The block- 
house walls, floor and ceiling were treated with urethane foam to pro- 
vide a semi-anechoic chamber within which acoustic intensity measure- 
ments were taken. Figure (4.2) details a cut-away front and top view 
of the chamber and test section. The acoustic treatment consisted of 
a base of 10cm foam, then a set of 10cm foam squares spread randomly 
and finally a set of 2cm foam squares also spread randomly. The 
effectiveness of the coating was determined by [14] using a horn 
driver and measuring the sound field decay with distance. The ane- 
choic treatment was effective between 200 Hz and 2000 Hz and only par- 
tially effective at 100 Hz. 
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Fiflure 4.1: Wind Tunnel FecUity 

4.2 Orifice/Louver Samples 

Nine separate orifice/louver samples were tested. All samples 
were made of plexiglass and screwed into the sample mounting window 
which was mounted flush into the side to the test section as shown in 
Figure (4.2). The separate samples tested were: 

(1) 0-1: a 4 1/2” (height) by 1 2/3” (width) rectangular ori- 
fice, 3/4” thick with straight sided walls 
the above made of 1/4" thick plexiglass 



(2) 0-2; 



I 

I 






I 

i 



40 . 




TOP View 
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(3) 

(^) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 



0-3: a 4 7/16" by 27/64" rectangular orifice, 1/4" thick 

with straight sided walls 

0-2A: 0-2 with 40 degree slanted walls, slanted away from the 
flow to reduce the downstream edge impact angle 
0-3A: 0-3 with 40 degree slanted walls 

L-1: 4 rectangular orifices, each 4 7/16" by 27/64”, separ- 

ated by 1/4", 1/4" thick 

L-2: 54 circular orifices, each 27/64" in diameter, separ- 

ated by 0.6328" between centers, 1/4" thick 
L-lA: L-1 with 40 degree slanted walls 
L-2A: L-2 with 40 degree slanted walls. 



The louver configurations were designed to maintain the same total 
open area as the single orifices, 0-1 and 0-2. The 40 degree slanted 
walls were slanted to provide a blunt trailing edge and subsequently a 
sharp leading edge. 



4.3 Test Section 



The test section was made from 3/4" finished plywood, 38cm by 38cm 
square in cross section and was located on the centerline of the cham- 
ber. The test section was sealed and the ends were further covered by 
2" of urethane foam. The background power levels, measured with the 
orifice covered, were sufficiently low compared to the power levels 
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with the orifice exposed (see Section 5.3) to ensure signal to noise 
ratios of at least 10 dB over the frequency range of Interest and at 
all operating speeds . 

4.4 Turbulent Intensity Measurements 

Turbulent Intensity is defined by the relationship: 



Turbulent fluctuations were measured using a Dantec Constant Tempera- 
ture Anemometer with a 5 micron diameter tungsten wire. The output of 
the anemometer, being the total voltage E, was inputted directly to 
the B6cK 2032 Spectrum Analyzer. Here the AC-coupled autospectrum, 
Gg(w), was measured. The mean-squared fluctuating velocity autospec- 
trum, G^(w) was deteirmined from: 



where H is the transfer function. See Appendix 2 for the details in 
determining H. Turbulent Intensity spectra were measured at all four 
flow speeds for orifice 0-2 and louvers L-1 and L-lA. The position of 
the U probe was : 




(4.1) 



G^(u>) - (|h|2) Ge(w) 



(4.2) 
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(1) two- tenths of an inch Into the boundary layer directly above 
the orifice at both the upstream and downstream edge (for the 
louver, the streamwise measurements were taken in the center 
of the furthermost upstream and downstream orifice, and 

(2) half-way into the throat of the orifice and at the same 
streamwise location as described above. 
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5 . EXPERIMENTAL RESULTS 



5.1 Introduction 

Figure 5.1 shows the geometry of the problem. Air flows at free 
speed Ueo over an orifice or louver which has been inserted flush into 
a long flat walled duct. The orifice was either rectangular, with 
streamwise dimension b and length L or circular with diameter d. 
The acoustic radiated power from the orifice/louver was measured at 
free stream velocities of 20, 30, 40 and 50 meters per second. The 
acoustic power was determined by sweeping the intensity probe over a 
control surface (a five sided boxed frame surrounding the orifice) 
then multiplying the resulting spatially and temporally averaged 
intensity by the total surface area of the control surface. Five ori- 
fices and four louvers were tested as described in Section 4.2. 







Figure 5.1: Sche^natic of orifice and probe geometry. Air flows at freestream velo- 

city, U, over an aperture of streamwise dimension, b. Sound intensity 
is measured by probe (microphone spacing, Ar) while continuously sweep- 
ing control surface. 
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To determine the contribution made solely by the orifice/louver 
source, two additional experiments were conducted. One experiment 
consisted of measuring the acoustic power radiated by the duct when 
the orifice was covered. The second experiment consisted of measuring 
the power flow through the orifice, in the absence of free stream 
flow, due to an artificially created acoustic field in the duct. This 
field was created using a speaker in the duct and was equal in magni- 
tude to that caused by the turbulent boundary layer at U - 50 m/s. 

This experiment determined the "open-window" contribution to the 
overall sound power. Finally, directivity measurements were taken for 
orifice 0-2 and louver L-2A at U - 25 m/s and the respective polar 
plots produced. 

5.2 Reactivity of the Sound Field and Sound Intensity Error Limits 

The reactivity of the sound field in the blockhouse (the closed 
acoustic space) was measured for various orifice and louver samples 
(see Figures 5.2 and 5.3). The Table below lists the measured dif- 
ference between the total sound pressure level, Lp and the total sound 
intensity level, for the orifice 02A and the louver L2A at a 

radial distance of 2 feet. 
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1 

1 

1 


Sarrple 


1 

1 

1 


Reactivity (Lp-L^ 


j U^(m/B) 


1 

1 

1 


1 

1 


02A 


1 

1 


0.8 


1 ^0 


1 

1 


1 

1 




1 


1.0 


1 50 


1 


1 

1 


L2A 


1 

1 


0.9 


1 50 


1 

1 


1 

1 




1 


1.0 


1 50 


1 

1 



Table 5.1: Measured Sound Field Reactivity 





Figure 5.2: Sound field reactivity for orifice 02A. Shown are the sound inten- 

sity and sound pressure spectral levels. The difference between the 
two is defined as the sound field reactivity. 
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Figure 5.3: Louver L2A sound intensity spectral levels versus sound pressure 

spectral levels (L^ vs. Lp). The difference between L^ and Lp is 
defined as the sound field reactivity. 



From Figure 2.4, the measured phase -mismatch for the sound inten- 
sity probe is: 



Frequency (Hz) 

63 

500 

1000 

2000 



(degrees) 

0.8 

0.2 

0.5 

0.9 



The maximum phase-mismatch between channels A and B for the Ana- 
lyzer is 0.3 degrees. With a spacer value of 50mm and the above sound 
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field reactivity, the phase reactivity nomograin of reference 6 gives 
the following value for the maximum phase of the sound field, between 
the two microphone positions. 



Frequency (Hz) 

63 

500 

1000 

2000 



(decreet) 

2 

16 

25 

60 



Therefore the limits of the phase -mismatch are: 



Frequency (Hz) 

63 

500 

1000 

2000 



<<^flrees) 

0.9 

15.5 

2A.2 

5S.8 



^degrees) 

3.1 

16.5 

25.8 

61.2 



The above total phase -mismatch values are the limits of error in 
the measured phase at the frequency listed. The corresponding limits 
in dB for the measured intensity from the nomogram are: 



Frequency (Hz) 



Error limits (dB) 



63 


♦ 2.5 


500 


♦ 0.5 


1000 


♦ 0.25 


2000 


♦ 0.125 



The corresponding statistical standard deviation, €, of the mea- 
sured spectral levels about the mean value is: 



I 



a 
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« - 4.34 / ySTa . (dB) (5.1) 

where B is the filter bandwidth and the total time record length. 

Substituting the appropriate values for B(2 Hz) and (150 sec) 

gives € - 0.25 dB. Thus the limits of error in measuring the sound 
intensity are approximately equal to the inherent statistical fluctua- 
tions (for frequencies above 100 Hz). 

5.3 Length and Frequency Scaling 

The three significant length scales for this experiment were: 

(1) the streamwlse orifice dimension, b 

(2) the turbulent boundary layer thickness, 6, above the orifice, 
and 

(3) the lateral correlation length, i. 

Table 5.2 provides the applicable values for b, S and i used in 
this experiment. 

The significant frequency scales were the non-dimensional strouhal 
numbers St(b) and St(i) defined by: 
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St(b/2) - f^b/Ueo (5.2) 

St(i) - cji/Uc - 





Orifice/I 
Louver | 


1 

b(cm) 1 
1 


1 

U(m/s) 1 
1 


1 (cm) i 


1 

<(cm) 1 

1. 


4/b 1 

for 0-2 1 


6/b 

for 0-3 


1 

1 

1 




0-1 1 


1 

i.23 1 

1 


1 

20 1 


2.4 1 


1 

4.19 1 
1 


0.99 1 


3.95 


1 

1 

1 




0-2,/ 1 


1 

^.23 1 


1 

30 1 


1.8 1 


1 

3.86 1 


0.91 1 


3.64 


1 

1 




0-2A 1 


1 

1 


1 

1 




1 

1 






1 

1 




0-3,/ 1 


1 

1.07 1 


1 

40 1 


0.9 1 


1 

3.65 1 


0.86 1 


3.44 


1 

1 




0-3A 1 


1 


1 

1 




1 

1 






1 

1 




L-h/ 1 


1 

1.07 1 


1 

50 1 


1.4 1 


1 

3.49 1 


0.83 1 


3.29 


1 

1 




L-IA 1 


1 


1 

1 




1 

1 






1 

1 




L-2,/ 1 


1 

1.07 1 


1 

50 1 


2.1 1 


1 

3.49 1 






1 

1 




L-2A 1 


1 


1 

1 




1 

1 






1 




Full 1 


1 

2.50 1 


1 

16.5| 




1 

4.70 1 


1.88| 




1 

1 




Scale 1 


1 

1 - 


1 

1 - 




1 

1 . 






1 

1 



Note: Full scale assones a louver positioned ISft from the leading edge of a 

flat surface in seawater at 58 F. The louver consists of Multiple rectan- 
gular shaped openings, each of transverse dimension, b^2.Scm. From Figure 
5.60 at 4/b » 1.88, the Strouhal # » 0.19. Therefore the full scale 
frequency, f^, is: 

f^ * (0.19)(16.5)/ (0.025) * 125 Hz 
Table 5.2 Experimental Length Scales 



where is the center frequency of the 1/3 octave band with the peak 
sound power level. Figure 5.4 is a graph of versus U^o for orifice 
0-1 at a constant b. The data points reflect several separate 
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experiments measuring the sound pressure level at radial distances of 
9, 18, 24 and 36 inches from the orifice. Figure 5.5 plots f^ versus 
for orifices 0-1, 0-2, 0-3 and louver L-1 and L-2. Clearly there 
exists a linear relationship (or nearly linear) between f^ and for 
a particular transverse dimension b. 
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Figure 5.4: Frequency/ velocity relationship for orifice 01. Shown are four sep- 

arate exp>eriinents for various radial distar>ces, r. 

Table 5.3 lists the non-dimensional Strouhal numbers for the above 
orifices/louver as well as the non-dimensional length scale, S/h for 
U - 20m/s. The strouhal number, St-fb/U, has been defined to equal 



{ 

I 

I 

I 

I 

I 



( 



C«nt«r fr«qu«ncy of p^ok 1 / J OB (Hi) 



52 - 



that recommended by Demetz and Farabee (5) to compare results. Figure 
5.6 is the graph of the data in Table 5.3 along with the results of 
Demetz and Farabee [5] for TBL flow over orifices enclosed by solid 
cavities . 




Figure 5.5: Frequer>cy/veloci ty relationships for orifices 01, 02, 03 and louvers 
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Orifice/ 

Louver 


l 

I St(b) . fb/U 

I Strouhal # 

I. 


1 

1 Non-dimensional 

1 length scale s/b 

1 


0-1 


I 

I 0.17 


1 

1 1.22 


0-2 


I 0.23 


1 1.22 


0*3 


I 0.12 


1 4.79 


L-1 


1 0.076 


1 4.79 


L-2 


1 0.048 

.1 


1 4.79 

J 



Table 5.3: Mon* dimensional Frequency and length Scales at U = 20 m/s. 




T9L Ihlcknesa/operiuro width, S/b 
□ OeMMz & Torobee + Conn 



Figure 5.6: Strouhal norter versus boundary layer/aperture width ratio 
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The decreasing trend in the Strouhal number with Increasing 6/b 
for the two cases agrees reasonably well. Having a knowledge of the 
boundary layer thickness, one can enter Figure 5.6 at the appropriate 
value for 6/b and determine the Strouhal number for orifice or louver 
flow noise. From the Strouhal number and the geometry of the apera- 
ture(s), one can predict f^. 

5.4 Structural Radiation Due to TBL Excitation 

To determine the contribution from structural radiation caused by 
the excitation from the TBL, an experiment was conducted at all free 
stream velocities with the orifice/louver covered. Identical sound 
power spectrums were recorded as described in Section 4.1. The 
resulting background power spectrums levels, Figures 5.7, 5.8, 5.9 and 
5.10, were 10 to 46 dB lower than with the orifice exposed. From 
these figures, the major contribution from the orifice/louver sources 
was in the range of 65 to 1000 Hz. Background noise was within 15 dB 
of the orifice/louver noise below 50 Hz above 1000 Hz. 
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Figure 5.8: Signal to noise for orifice 02 at U * 40 m/s. 



I 





POWER. D0 RE lE-12 WATT POWER. OB RE lE-12 WATT 





0. 500. 1000. 1500. 2000. 

FREQUENCY (HZ) 



Figure 5.10 Signal to noise for louver 11 at U = 50 m/s. 
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5.5 "Open Window" Contribution 

Previous work by Spinka [15] determined the acoustic pressure 
field within the test section at various free stream velocities. To 
determine the contribution from the open window affect, a similar 
acoustic pressure field was generated within the duct using a speaker 
driven a B&K random noise generator, type 1402, and a B&K band pass 
filter set, type 1612. The resulting power which radiated through the 
orifice was measured as described in section 2.5 and compared to the 
power radiated by the orifice source. Table 5.4 lists the resulting 
power levels for both cases at U - 50m/s. The conclusion was that the 
open window affect results in levels which were typically 23dB lower 
than the levels due to the orifice source. 



1/3 oe center 
frequency 



Open window 
power in 1/3 06 (2) 



1 

400 1 


53.8 


1 

1 


77.6 


1 

1 23.8 


500 1 


48.2 


1 


72.7 


1 24.5 


630 1 


45.9 


1 


69.2 


1 23.3 


800 1 


43.7 


1 


67.9 


1 24.2 


1000 1 


46.2 


1 


66.4 


1 20.2 


1250 1 

1. 


41.0 


1 

1 


64.5 


1 23.5 

1 



Orifice 0*1 
power in 1/3 06 (3) 



(3-2) 



Table 5.4: "Open window" power levels versus orifice 0 -1 power levels at U = 50m/s 
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5.6 Orifice Source Power Spectra and Total Power 

The measured sound power spectra versus the free stream velocity, 
U-20, 30,40 and 50 m/s, for 0-1, 0-2, 0-3, L-1 and L-2 are shown in 
Figures 5-11, 5-lla, 5-llb, 5-llc, 5-lld. The total sound power radi- 
ated by each orifice/louver source at the above speeds is given in 
Table 5.5. 




Figure 5.11 Orifice 01 sooxl power spectral level versus free stream velocity 
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Figure S. 11(a): Orifice 02 sound power spectral levels versus freestrcam velocity 
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Figure 5.11(b): Orifice 03 sound power spectral levels versus freestream velocity 
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Figure 5.11(c): Louver Li sound power spectra levels versus freestream velocity 
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Figure 5.11(d): Louver L2 sound power spectral levels versus freestream velocity 
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1 

1 

1 

1 


Orifice/ 

Louver 


Free Stream 
Velocity U(m/s) 


I Total Sound Power 
I <dB re 1 E*12 watts) 




1 

1 


0-1 


20 


1 63.3 




1 




30 


1 73.3 




1 




40 


1 79.6 




1 




50 


1 84.0 




1 

1 


0-2 


20 


1 62.7 




1 




30 


1 73.2 




1 




40 


1 78.4 




1 




50 


1 83.0 




1 

1 


0-3 


20 


1 62.9 




1 




30 


1 72.2 




1 




40 


1 76.3 




1 




50 


1 79.7 




1 

1 


L-1 


20 


I 74.5 




1 




30 


1 82.5 




1 




40 


1 89.2 




1 




50 


1 94.5 




1 

1 


L-2 


20 


1 64.9 




1 




30 


1 73.6 




1 




40 


1 80.1 




1 




50 


I 85.3 




1 

1 


0-2A 


20 


I 62.0 




1 




30 


I 71.8 




1 




40 


1 77.2 




1 




50 


1 82.2 




1 

1 


L-1A 


20 


1 58.1 




1 




30 


1 68.4 




1 




40 


1 73.3 




1 




50 


1 77.7 




1 

1 


L-2A 


20 


1 55.3 




1 




30 


1 65.0 




1 




40 


1 69.7 




1 

1 




50 


1 76.3 





Table 5.5: Orifice/louver total sound power levels versus freestream velocity.U, 
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Several observations are made from reviewing the above data. 

5.6.1 The effect of louvering on sound power 

(a) louver, L-1, radiated the most acoustic power, averaging lOdB 
more at all speeds than the single orifice though the total 
open areas were the same (see Figures 5.12 and 5.12a). 

(b) louver, LI, showed strong resonant peaks and had the narrow- 
est half -power bandwidth of all the samples tested. 

(c) louver, LI, radiated significantly more power than louver L2 
in the lower frequency bands but approximately equal power in 
the middle to high frequency bands (see Figures 5.13 and 
5.14) 




Figure 5.12 The effect of louvering on power, L1 versus 02 at U=30iVs 
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Figure 5.12«: The effect of louvering oo sound power, LI versus 02 at U=50m/i 
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Figure 5.13: louver geometry effect on sound power, 11 versus at L2 at U = 30m/s 
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(d) louver, L2 , had the broadest spectra of all the samples 
tested with equal open area. 

(e) louver, L2's total radiated power was within 2 dB of orifice 
02 's power in all speed ranges but the spectral shape was 
significantly broader. As shown in Figures 5.15 and 5.16, 
louver L2 radiated significantly more power than 02 in the 
middle to higher frequency bands (10 dB more on average) but 
slightly less (approximately 1-2 dB) power in the lower 
frequency bands. 
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Figure 5.15: Effect of louvering on sound power, L2 versus 02 at U = 30fn/t 




Figure 5.16: Effect of louvering on sound power, L2 versus 02 at U = 50 m/s 
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Thus the effect of louverlng, simply by dividing the total open area 
into more numerous, smaller areas does not by itself reduce the total 
radiated flow noise and will most probably result in increased radi- 
ated sound power in the higher frequency bands. 

5.6.2 Orifice throat thickness effect on radiated power 

Figure 5.17 shows the orifice throat depth effect on sound power 
by comparing orifice 01 to orifice 02. Orifice 02 radiated more power 
(approximately 3dB) in the higher frequency bands but maintained 
approximately 1-2 dB less power in the more dominant lower frequency 
bands . 




Figure 5.17; Orifice depth effect on saxd power, 01 versus 02 at U = 30 and 50 m/s 
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Thus the throat thickness was more predominant In effecting the 
spectral shape than reducing the total radiated power. The effect of 
reducing the orifice throat thickness appears to broaden the spectral 
shape while maintaining the same total radiated sound power. 

5.6.3 Streamwise dimension, b, effect on sound power 

The spectral shape of orifice 03 had the widest half-power band- 
width. The effect of reducing the streamwise dimension, b, was to 
broaden considerably the power spectra while reducing only moderately 
the total sound power. From section 3, the controlling shape factor 
was (l-cos(Kb)). Thus the predicted result from reducing b by four 
fold would be a subsequent increase by four fold in the power spectra 
bandwidth. Reviewing again the power spectra for orifices 02 and 03 
(Figures 5.11a and 5.11b), we see that indeed the power bandwidth for 
03 is approximately four times that of 02 (defining the bandwidth to 
be between the one-quarter peak spectral power points). 

5.7 Speed/Power Laws 

The speed/pover law is of the form: 

P - C 

where C is a value which is independent of speed. The data in Table 
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5.5 was utilized to determine the total radiated power dependence on 
speed. The equations used to determine the value of the speed expo- 
nent, n, were: 

P3/P2 - (U3AI2)'' 

10 log (P3/Pref) - 10 log(P2/Pref> “ n 10 log(U3/U2> 

Lp3* Lp^- n log(U3/U2> 

n - (Lp^- Lp^/(10 log (U 3 /U 2 )) 

where the subscripts 3 and 2 correspond to 30 and 20 m/s. Similar 
equations hold for U 4 /U 2 , U 5 /U 2 , U 4 /U 3 , U 5 /U 3 and U 5 /U 4 . Table 5.6 
lists the calculated value for n for each sample tested. The values 
for n varied from 3.3 to 6.8 with a mean value of 4.9 and a standard 
deviation of 0.3. Excluding the extremely low values of n recorded 
for 03 at the higher free stream velocities, the mean value of n was 
5.0 with a standard deviation of 0.2. 

The speed/power relationship for orifices 01, 02 and louvers LI 
and L2 is plotted in Figure 5.18. Also shown is the least squares fit 
to the data points satisfying the equation: 

(Lpi - LPj) - n 10 log(Ui/Uj) (5.2) 



From the fit, the best value for n is 4.97. 
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Orifice/Louver 

0-1 
0-2 
0*3 
L-1 
L-2 
0-2A 
L-U 
L-2A 
average 
std. dev. 



U3/U2 U4/U2 

5.7 5.4 

6.0 5.2 

5.3 4.5 

4.6 4.9 

4.9 5.2 

5.6 5.1 

5.9 5.1 

5.5 4.8 

5.4 5.0 

0.5 0.3 



U 5 /U 2 


U 4 /U 3 


5.0 


5.2 


4.2 


5.1 


3.3 


4.2 


5.4 


5.0 


5.5 


5.1 


4.3 


5.1 


3.9 


4.9 


3.8 


5.3 


4.5 


5.0 


0.6 


0.3 



U5/U3 U5AJ4 

4.8 4.5 

4.4 4.7 

3.4 3.5 

5.4 5.5 

5.3 5.0 

4.7 5.2 

4.2 4.5 

5.1 6.8 

4.7 5.0 

0.6 0.9 



overall average 4.9 overall standard deviation 0.3 



5.6 Speed/Power relationship (P r o/'), values for the speed exponent, n 
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Figure 5.18: Speed/Power relationship for four samples (P=CU*^) with least squares fit 
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Thus the concluding dependence is: 

Total radiated power increases with the 5th power of the free 
stream velocity 

5.8 Directivity of orifice 01 and louver L2A 

Directivity measurements from orifice 01 and louver L2A were 
taken at U - 25m/s. Figure 5.19 details the coordinate system used 
for the measurements. The measurements indicated the source was 
nearly omni-directional. Power levels in the XZ (azimuth) plane dev- 
iated by 0.9 dB from the maximum level. In the YZ (polar) plane, power 
levels deviated by 2.6 dB from the maximum level. Off the major axis, 
power levels deviate by 3.6 dB from the maximum level directly above 
the orifice source. Table 5.7 lists the measured sound intensity in 
dB re IE- 12 watt/m^ at a radial distance of r - 27 inches and the 
various polar and azimuth angles. Figure 5.20 is the azimuth direc- 
tivity plot for orifice 01 and louver L2A and clearly supports the 
conclusion of a nearly omni-directional source. 
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Figure 5.19 Coordinate systefn for directivity measurefnentt. Reference is the 
aperture center. Polar angle, $, ineasured from 2 axis, azimuthal 
angle, measured from -x axis. Measurements described by spheri 
cal coordinates, r,$,4. 



Total Intensity 



lar angle 


Azimuthal ar>gle 


dB re 1E-12^ 




Orifice 0-1 




0 


90 


57.3 


0 


45 


56.9 


0 


15 


56.4 


0 


135 


57.1 


0 


165 


56.8 


0 


105 


57.2 


0 


75 


57.0 


•85 


90 


54.9 


•45 


90 


56.5 


0 


90 


57.5 


45 


90 


57.0 


85 


90 


55.9 


•45 


30 


54.7 


•45 


25 


53.9 


-45 


35 


55.0 



Table 5.7 Directivity measurement for orifice 0-1 and louver L2A 
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Total Intensity 






Polar angle 


Azimuthal angle 


d8 re 1E*12watt/m^ 








Louver L2A 








0 


160 


66.2 






0 


145 


66.1 






0 


130 


66.0 






0 


115 


66.0 






0 


90 


65.8 






0 


75 


65.7 






0 


60 


65.7 






0 


45 


65.8 






0 


30 


65.2 






0 


15 


65.2 






•90 


90 


63.4 






•75 


90 


63.9 






•60 


90 


65.0 






•45 


90 


65.4 






•30 


90 


65.9 






•15 


90 


66.1 






0 


90 


66.2 






15 


90 


66.4 






30 


90 


66.5 






45 


90 


66.6 






60 


90 


66.4 






75 


90 


66.3 






90 


90 


65.6 






•45 


45 


65.3 






•45 


30 


64.2 






•45 


90 


65.3 






•45 
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66.0 





Table 5.7 Directivity measurefnent for orifice 0*1 and louver L2A (Continued) 



73 




Figure 5.20 Orifice 01 directivity in the XZ (azimuth) plane and louver L2A 
directivity in the Y2 (polar) plane 
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5.9 Edge Angle Effect on Radiated Sound Power 
5.9.1 Introduction 

Earlier it was shown that Indeed the orifice source provided both 
the spectral character and the predominant amount of the total power. 
The strength of this source is due in part to the eddy interaction 
with the downstream edge of the orifice. To determine the relative 
importance of this edge interaction compared to the other significant 
parameters identified in Section 3, three additional orlf ice/louver 
samples were tested. They were orifice 0-2A and louvers L-lA and 
L-2A. The details of these samples are described in Section 4. The 
only difference between these samples and the previously tested 
samples was the leading and trailing edges were slanted away from the 
flow by 40 degrees. This provided both a blunt trailing edge where 
the eddies strike and a sharpened leading edge (essentially reducing 
the upstream aperture to 1/32", the radius of the edge). 

5.8.2 Angled edge orlf Ice/louver power spectra 

The measured sound power spectra for the above samples are shown 
in Figures 5.21, 5.22 and 5.23. The total power radiated by each 
sample at the various freestream velocities is given in Table 5.5. 
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Figure 5,21 : Orifice 02A power spectra level versus freestreao velocity 




Figure 5.22: Louver L1A power spectral levels versus freestream velocity. 
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Figure S.23: Louver L2A power spectral levels versus freestream velocity. 

Several observations are worth noting: 

(a) the effect of angling the edges is more predominant for the 
louvers than the single orifice. 

(b) the effect on louver L-IA is a 14-16 dB reduction in total 
radiated sound power. The large resonant peaks which occurred 
in the spectra for L-1 are not present in the spectra for 
L-lA. Figures 5.24 and 5.25 compare the two louvers at U=20 
and 50 m/s. Not only are the resonant peaks eliminated by 
angling the edges but also the radiated power is reduced by 
5-7 dB in the higher frequency bands. 
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Figure 5.24: Edge angle effect on sound power spectral level, LI versus LIA 



at U s 20 m/s 




Figure 5.25: Edge angle effect on power spectral level, Li versus L1A at 
U = 50 m/s 
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(c) Figure 5.26 compares the power spectral levels of L-2 and 
L-2A. The result for this louver is a 8.6 to 10.4 dB reduc- 
tion in sound power equally divided across the frequency 
bands . 

(d) for the single orifice, 0-2A, the effect of angling the edges 
resulted in only a 0.7 to 1.4 dB reduction in sound power. 
Figures 5.27 and 5.28 compare the power spectra for orifice 
0-2 and 0-2A at U-20 and 50 m/s. Here again the effect of 
angling the edges resulted in a reduction in sound power in 
the higher frequency bands but in this case, essentially no 
change in the lower frequency bands . 




Figure 5.26: Edge angle effect on sound power spectral level, 12 versus L2A at 
U = 50 m/s. 
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Figure 5.27: Edge angle effect on sound power spectral level, 

02 versus 02A at U ^ 50 m/s. 
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Figure 5.28: Edge angle effect on sound power spectral level, 

02 versus 02A at U = 20 m/s. 
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5.9.3 Combined effect of louvering and angling edges 

Taken together, the effect of louvering and angling edges results 
in a reduction of 3.4 to 7.5 dB in radiated sound compared to the 
single angled orifice (both of equal total open area) . This reduction 
most probably is caused by two factors: 

(a) the louvers shed eddies which after convecting downstream to 
the trailing edge are the same order in size as the edge 
thickness. This results in a more beneficial interaction when 
the eddies impact the angled edge. The orifice, however, with 
a four times larger aperture, sheds eddies which are much 
larger in size upon impact with the trailing edge and are 
therefore less affected by the thin angled edge. 

(b) louvering, by reducing the transverse dimension, b, interferes 
with the longitudinal correlation of the shed eddies. This 
interference will be more predominant in the lower frequencies 
where the longitudinal correlation length is larger than the 
aperture transverse dimension, b. As shown in Figures 5.29, 
5.30, 5.31 and 5.32, the reduction in sound power is predom- 
inantly in the lower frequency bands. There is, however, more 
radiated power in the angled louver in the higher bands by 
approximately 1-3 dB. 
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Figure 5.29: Effect of Louvering on socnd power spectra level, L2A 

versus 02A at U ^ 30 m/s 




Figure 5.30: Effect of louvering on sound power spectral level, L1A 

versus 02A at U = 30 m/s 
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Figure 5.31: Effect of louvering on sound power spectral level, l2A versus 02A 

at U = 50 m/s 




Figure 5.32: Effect of louvering on sound power spectral level, LlA 

versus 02A at U = 50 m/s. 
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Thus one can conclude that effective louvering does reduce the 
total radiated power by reducing single mode resonance and by inter- 
fering with the lower frequency, longitudinal eddy correlation. This 
results in effectively broadening the power spectrum. 

5.9.4 Louver geometry effect on sound power 

As with the straight walled louvers, the nniltiple circular aper- 
ture louver was quieter than the multiple rectangular aperture louver 
by 1.4 to 3.6 dB in total radiated power. As shown in Figures 5.33 
and 5.34, louver L-2A had approximately the same spectral shape as 
louver L-lA but with considerably less power in the low frequency 
bands. Again this reduction can be attributed to the additional 
interference of the lateral correlation due to the reduced dimension 
in the aperture's height. To test this conclusion, one can calculate 
the frequency where the lateral correlation length, i, equals the 
aperture diameter, d for louver L-2A and at freestream velocity of 50 
m/s. From Figure 3.4, at B(K(r3))-0.4, the value of Kr 3 equals 1.2. 
Setting r 3 equal to 0.01 m (27/64") and U-50 m/s gives a limiting 
frequency of 450 Hz. Thus one would predict diminishing reductions in 
sound power for frequencies above 450 Hz. Reference to Figure 5.34 
shows that, indeed, at frequencies above 450 Hz the effect of lateral 
correlation interference on reducing sound power is diminishing. 



60 . 



84 



c\j 

I 



UJ 

q: 

Q 



q: 

UJ 

:» 

o 

a. 



< 



I 

UJ 



UJ 

cc 

CD 

a 



cr 

UJ 

o 

a. 




0 - 500 . 1000 . 1500 . 2000 . 

FREQUENCY (HZ) 

Figure 5.33: Louver geometry effect oo sound power spectral level, 

LIA versus L2A at U > 30 m/s. 




Figure 5.34: Louver geometry effect on sound power spectral level, 

L1A versus L2A at U = 50 m/s. 
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Figure 5.35: Turbulent intensity spectral levels for louver LI versus freestream 
velocity. 



5.10 Describing the Velocity Field Surrounding the Samples 

In Section 3, the radiated sound power was shown to be directly 
related to the normal component of the turbulent intensity. Several 
measurements of the turbulent intensity spectrum were taken as 
described in Section 4.4, The assumption was made that the turbulent 
intensity in the streamwise and normal component directions were 
equal. Thus 3 dB was subtracted from the measured intensity spectral 
levels to give only the normal component of the turbulent fluctua- 
tions. The following observations are made: 
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(a) the turbulent Intensity measured 2 tenths of an inch Into the 
boundary layer was relatively unaffected by the probe loca- 
tion in the transverse direction and the various orifice 
geometries. Figure 5.35 is the measured normal component of 
the turbulent intensity for louver L-1, with the probe at the 
center of the upstream aperture. These spectral levels were 
used in the following section to collapse the sound power 
spectra. 

(b) there were significant reductions in turbulent intensity 

(approximately 3 dB ) between louvers L-1 and L-lA when meas- 
ured at the leading edge of the leading aperture. This reduc- 
tion in intensity supports the measured reduction in total 
sound power and also supports the conclusions that the sound 
power is related to the driving aerodynamic field. 

Despite a fairly rigorous attempt to maintain the block-house air 
tight, a finite and measurable pressure differential existed between 
the blockhouse and the duct static pressure. This pressure differen- 
tial resulted in a small amplitude, mean flow entering the ori- 
fice/louver from the closed acoustic space. The magnitude of this 
transverse mean flow was greatest in the louver samples and measured 
on average, fourteen percent of the freestream velocity. The exact 
contribution of this mean flow to the radiated sound power in not 
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known but it Is estimated to be of minimal Impact since the turbulent 
Intensity of the entering flow measured several orders of magnitude 
less than the normal component of the freestream turbulence. 

5.11 Verifying the Scaling Laws 

In Section 3, an equivalent source model was developed for rectan- 
gular shaped aperture orifices and louvers. To ascertain the validity 
of this model, the sound power spectra for orifices 02, 03 and the 
louver LI were collapsed by normalizing the spectral levels by the 
magnitude of the equivalent source model. The resulting collapsed 
spectra are shown in Figures 5-36, 5-37 and 5-38. It is evident by 
the collapse that the equivalent source model adequately describes the 
controlling parameters for orifice and louver flow noise. The predic- 
tion of the actual magnitude of the sound power is fairly accurate for 
the samples with the smaller apertures but is too high by about a fac- 
tor of 100 for the larger aperture openings (01 and 02). The probable 
reasons for this are (1) overestimating the effective aperture width, 
b, (2) the estimation for the lateral correlation lengths and (2) the 
assumption of equal components of turbulent intensity. As mentioned 
in Section 3, for orifice 0-2, the longitudinal correlation amplitude 
function, A(K(ej^) does attenuate to a value of approximately 0.6 for 
the slower speeds indicating that some eddy decay occurred in the 
aperture. For this case, one could assume a smaller transverse 
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length, b' , over which the eddy decay was raintmal and use this length 
scale for evaluating the source strength. Defining b' at the point 
where the magnitude of the amplitude function equaled 0.95 results in 
an effective aperture length of 0.5 inches versus the actual dimension 
of b-1.67 inches (for orifice 0-2 at U-20m/s and f-250 Hz). This 
results in a reduction of 4.2 dB in the predicted spectral level and 
could account for some of the error in predicting the larger aperture 
radiated sound power. 




Figure 5.36: Scaling law verification for orifice 02 
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Figure 5.38: Scaling law verification for louver Li 
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6, CONCLUSIONS 

The acoustic source created by turbulent flow over orifices and 
louvers was Investigated by measuring the radiated sound intensity, 
utilizing an intensity probe, at various wind speeds in a semi- 
anechoic wind tunnel. From the intensity spectra, the radiated power 
and the directivity of the source was determined. The goal of this 
research was to describe the source and identify the controlling par- 
ameters which affect its spectral shape, radiated power and directi- 
vity. The orifice and louvers investigated have geometries similar to 
those used in industry in both air and water environments. The 
effects of various leading/trailing edge geometries were also investi- 
gated. 

The source of noise has been modelled on the basis of the theoret- 
ical analysis by Ffowcs Williams, Nelson and Corcos. The equivalent 
source of noise was modelled as two equal but opposite monopole 
sources located in the upper and lower aperture of each orifice. The 
strength of each source is determined by the mass flow driven through 
the apertures by the normal component of the turbulent fluctuations. 
The resulting radiated sound field, both into a closed acoustic space 
and radiating away from the orifice, is dependent on the acoustical 
properties of the louver. 
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To ascertain the validity of the source model, two orifice and one 
louver power spectra were collapsed by normalizing the spectral levels 
by the magnitude of the equivalent source model. It is evident by the 
collapse that the equivalent source model adequately describes the 
controlling parameters for orifice and louver flow noise. The con- 
trolling parameters for the geometries tested were: 

(1) Strouhal number based on the lateral correlation length, i, 
and aperture transverse dimension, b 

(2) the normal component of the Turbulent Intensity, 

(3) the fifth power of the freestream velocity and 

(4) the shape factor (1 -cos (Kb)) where K is the convective 
wavenumber, K - c*?/Uc, and Uc is the eddy convection velocity. 

As a result of this research, several important conclusions are 
made : 

(a) sound intensity measurements , using an intensity probe, can 
successfully describe an acoustic source, even though 
measurements are taken within the source's near field. 

(b) the errors associated with measuring intensity, using an 
intensity probe, can be as small as the inherent statistical 
deviation of the random fluctuations . 
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(c) the equivalent source model developed to describe the behav- 
ior of orifice and louver flow noise successfully describes 
the controlling parameters and provides valuable insight for 
predicting the actual magnitudes of the radiated sound spec- 
trum. 

(d) the 1/3 octave band containing the most acoustic power from 
orifice and louver flow noise can be predicted from the non- 
dimensional Strouhal number based on the aperture's trans- 
verse dimension and the non-dimensional length scale, S/\> 
where S is the turbulent boundary layer thickness. 

(e) for aerodynamic flows, the orifice/louver source provides 
both the dominant amount of total sound power and the spec- 
tral character of the radiated noise. Contributions due to 
structural radiation from the ducting and ducting "open win- 
dow" effects are not significant. 

(f) the effect of louvering, simply by dividing the total open 
area into more numerous, smaller areas does not by itself 
reduce the total radiated flow noise and will most probably 
result in increased radiated sound power in the higher 
frequency bands. 

(g) the effect of reducing the aperture throat thickness is to 
broaden the spectral shape of the radiated power spectrum 
while maintaining approximately the same total radiated sound 



power. 
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(h) the effect of reducing the aperture's transverse dimension 
can be adequately described by the controlling shape factor 
in the equivalent source model, (l-cos(Kb)). 

(i) directivity measurements show a nearly omni-directional 
source in the on-axls planes with off-axis power levels 
reduced by 2-3 dB. 

(j) taken together, the effect of louverlng and angling the lead- 
ing and trailing edges resulted in reductions to sound power 
of 3.4 -7.5 dB compared to the equivalent area, single aper- 
ture, angled orifice. 

(k) the measured reductions in sound power due to effective lou- 
vering can be explained by (1) reductions in the normal com- 
ponent of the turbulent intensity due to the eddies impacting 
a blunt trailing edge and (2) destructive interference with 
the eddy longitudinal and lateral correlations. 

(l) the louvers with numerous circular apertures and angled edges 
radiated the least acoustic power with a very broad, flat 



power spectrum. 



0 
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APPENDIX 1: MATCHED ASYMPTOTIC EXPANSION SOLUTION TO A RADIATING 
ORIFICE IN AN INFINITE RIGID WALL 

The two dimensional problem to be solved is the far field radiated 
pressure at xj^, created by a driving aerodynamic field at y^, within a 
finite size aperture in an infinite, rigid wall. Using the principle 
of reciprocity, the above field can be determined by first solving for 
the sound field at yj^ created by a distant monopole at xj^. Then, the 
desired solution is obtained by simply replacing the position vectors 
of the far field monopole and the aperture source. The solution to 
this problem will be analyzed using the principles of Matched Asymp- 
totic Expansions (MAE). 

The coordinate system used in the analysis is shown below, 




with boundary conditions, 3p/9y2 “ 0 for |yxl - ^ 72 



The governing equations for a plane wave of single frequency, w. 



at normal incidence are: 



I 



! 

I 



I 
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V^p 4- k^p - 0 (Al . 1) 

i(ky^-wt) 

p - Ae (Al.2) 

where k is the acoustic wavenumber, k - <*>/co- Also define the 
non-dimensional terms, r - wt, - kyj^, y 2 - ky 2 , p - p/A and e - kb 
for the parameters in the outer problem. 

Substituting the above outer parameters into equations Al . 1 and 
Al.2 gives the governing equations and boundary condition for the 
outer problem 

v2p + p - 0 (A1.3) 

i(y2-0 

p - e (Al.4) 

To lowest order, €->0, and the orifice closes. Thus the outer 
problem becomes 

p - ^iPl ^2P2 *** (A1.5) 

with 6]^ - 1. From the boundary condition, we can determine p]^ to be 

iy2 -iy2 

+ e 



Pi - Pi Pr - ® 



, y2 < 0 



(A1.6) 
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suppressing the e'^^ factor, where p^^, pj. are the incidence and 
reflected waves, respectively. Note that the boundary condition is 
satisfied 



3pi ly2 -iy2 

- i (e -e ) - 0 for y 2 - 0 (Al.7) 

3y2 

Choose also the condition that p^ ■ 0 for y2 > 0* To open the 
orifice, choose the inner variables yj^ - y\/€ and y 2 - y2/^ • 

In Inner variables, 



p - exp(l(ey2-f) ] + exp[ -i(€y2-f ) ] 



- exp(-if ){exp(ey2> + exp(-€y2>) 



2~2 2“2 

- _ _ ^y 

- exp(-lf ) {l+ey2 - ^ * •••+! - ey2 - ^ +•••} 

2 2 

2-2 

- ^ y _ 

- 2e"^^ {1- 2 + 0(^^)) as € -♦ 0»y2 fixed (y 2 0) (A1.8) 

2 

The governing equation with the inner parameters substituted into 
equations A1.3 becomes 



v2p + <2p _ 0 (Al.9) 

with boundary conditions 3p/3y2 " 0 for y 2 ” 0 and lyj^l i 1. Let 
p be defined by 



P ~ Ao(<)pi + Ai(<)p2 + 



• • • 



(Al.lO) 



I 



[ 



i 



99 



for the inner expansion. The resulting inner problem coordinate sys 
tern is 







Since V^pj^ - 0, solve the inner problem by conformal mapping, pj^ - 
A Rew. Using a solution of the form, iz - cosh w, (suppressing the 
factor) we can determine 

1 1-3 

w - in 2 2 + i^r/2 + +••• as |z|-^ (Al.ll) 

J. 

2x2z 2x4x4z 

where z - re^^. 



Since z - z/€ , 



w - in 2z - ine + iir/2 + for € -♦ 0, z fixed (z-^c) (Al.12) 



Therefore, utilizing equation A. 11, we can determine Aq(€)p]^ to be 



AAQ{in 2r - ine + i7r/2 + ®(€^)} 



(A1.13) 



I 

I 



\ 

i 



I 

( 

I 



( 
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Now matching the leading terms in the outer and inner solutions, 
equations Al.8 and Al.13, leads to Aq - + 1/in € and A - - 2. Then 
S2(«) “ 1/in € in the outer expansion. P 2 satisfies the Helmholtz 
equation with k - 1 and must behave like In r for r 0, € -► 0 and r 
fixed. Clearly 



( 1 ) 

P 2 - C Hq (r) 



(A1.14) 



satisfies this requirement and also satisfies the rigid wall condi- 
tion. As r 0, 

( 1 ) 21 21 
Hq (r) - — inr + l + — (y - in2) + *** 
n n 



where y is defined as Euler's constant. Thus to match with the inner 
solution, A - - 2 - C • 2i/ir, C must be equal to ix. Therefore the 
final outer expansion becomes 



P 



.-ii' 



y2 *y2 iirHp)(r) 

(e + e + + • • • ) 

ine 



(A1.15) 



where the first term in brackets is the incident wave, the second term 
is the reflected wave and the lowest order orifice effect is the last 
term which describes a monopole type radiation. Substituting for the 
non-dimensional terms and applying the principle of reciprocity, one 
obtains the radiated sound from a 2 dimensional, single slot, line 
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source 



- kx2 -kx2 iwH^^^(r) -i«t 

P - A(e + e + )e (Al.16) 

inkb 



(1) - I i(r-jr/4) 

As r -«o, (r) -♦ j2/jrr e so the for field radiated 



sound for the 2D, single slot, line source becomes 



^ kx 2 -kx 2 I =-i i(r*jr/4) -iwt 

P - A(e + e + ijrj2/jrr e )e (A1.17) 

r-Ko 



Equation Al.17 provides the result that the far field radiated 
pressure for a 2D, single slot source to lowest order is composed of: 

(1) the incident and rigid wall reflection 

(2) 2D line source having a 1/r dependence (cylindrical spread- 
ing) . 

Thus one can infer that it is plausible for the 3 dimensional, 
single orifice source to have a reflection coefficient, |r|^ = 1 for 



values of b/A « 1. 
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APPENDIX 2: MEASURING TURBULENT INTENSITY 

In this appendix, the transfer function to relate the mean and 
fluctuating component of the hot wire anemometer output voltage to the 
mean and fluctuating component of the particle velocity is developed. 

The basic relationship between output voltage, E, and particle 
velocity, U, is defined by King's Law 

E - A + BU^ (A2.1) 

Now define the instantaneous voltage and particle velocity by 

E - Eq + e'. |e' I « |Eo| (a2.2) 

U - Uq + u' , |u' I « |Uol (A2.3) 

where Eq is the mean of E, Uq is the mean of U, e' and u' are the 
fluctuating components whose expectations are zero. Substituting 
equations A2.2 and A2 . 3 into equation A2.1, expanding and dropping the 
higher order terms of u' and e' gives 

e 2 + 2Eoe* - A + BU^ (1 + (nu'/Uo)) (A2.A) 

Taking the expectation of equation A2.4 gives 
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e2 - a + BU^ (A2.5) 

0 0 

Subtracting equation A2.5 from equation A2.4 gives the desired 
relationship between the fluctuating components of the anemometer 
output voltage, e', and the particle velocity, u' , 

1-n 

u' - ([2 EoUq ]/Bn)e' (A2,6) 

Thus, defining the term in parentheses as H, the resulting rela- 
tionship between the spectral densities of the fluctuating velocity, 
S^(f), and the output voltage, $£(f), is 



S^(f) - |h|2 SE(f) 



(A2.7) 
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APPENDIX 3: MATHEMATICAL STEPS IN DETERMINING SQ(yi^,cj) 

From equation 3.7, the cross spectral density is assumed to be of 
the form 

- S^(cj)A(Kci)B(Ke 3 ) • exp (iKci) (A3.1) 

As described in Section 3.2, the normalized longitudinal cross- 
spectra amplitude function, A(Kc]^), is identically equal to one for 
all samples except 01 and 02 at the lower freestream velocities, where 
at a minimum it equals 0.6. So assuming no decay in the eddy struc- 
ture in the yj^ direction (l.e., A(K€j^) equals one) and to be 

uniform over the aperture surface, the source strength spectral den- 
sity becomes 

Sq(yK.“). “ P^w^S^(yij,w)JJ'B(K< 3 )exp(iK«i)tf(«j^)dA(«i^) (A3. 2) 

00 

where €3 are the position coordinates as shown below. 
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Substituting for have 



SqCyic.w) 



pV s^(y^.w) J / B(K«3) 

CO CO 



ik€^ 

e (L- I < 3 ! ) (b- I «i I ) 



(A3. 2a) 



• dA («k) • 



Since the lateral correlation i, is typically one- tenth or less 
the length of the aperture, the Integral in the €3 direction can be 
evaluated, independent of the actual limits of the integrand, as 



I B(K€3)(L-|£3|)d«3 - Li - 
00 

where 



B(K«3)d£3 - i 



(A3. 3) 



(A3. 3a) 



Evaluating the Integral over the €j directions gives: 



‘ 1K£, 2 

e (b-|«]^|)d£i (l-cos(Kb)) 

J 

GO 



(A3. 4) 



Therefore the final result, neglecting the i^ term is: 



SQ(yk.w) - P^w^Sv(yk.«)Li(— ) (l-cos(Kb)) 






(A3. 5) 



I 



I 
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The remaining term needed to fully describe the magnitude of the 
source strength spectral density is i. The lateral correlation length, 
i, is determined at the point where the amplitude function, 

B(K€ 3 > equals 0.4 and the frequency equals the center frequency, f ^ , 
as defined in Section 5.3. 
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